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We introduce another powerful method of solving PDEs. First, we need to consider 
some preliminary definitions and ideas. 

1 Preliminary ideas and motivation 
1.1 The delta function 

Ref: Guenther & Lee §10.5, Myint-U & Debnath §10.1 

Definition [Delta Function] The 5-function is defined by the following three 
properties, 



where / is continuous at x — a. The last is called the sifting property of the (^-function. 

To make proofs with the 5-function more rigorous, we consider a 5-sequence, that 
is, a sequence of functions that converge to the (^-function, at least in a pointwise 
sense. Consider the sequence 
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Definition [2D Delta Function] The 2D 5-function is defined by the following 

three properties, 

0, (x,y)^0, 
oo, (x,y) = 0, 



5 (x,y) 





5 (x, y) dA = 1, 
/ (x, y) S (x - a, y - b) dA = f (a, b) . 



1.2 Green's identities 

Ref: Guenther & Lee §8.3 

Recall that we derived the identity 




f (GV-F + F- VG)dA= [ (GF)-hdS 

D Jc 



(1) 

for any scalar function G and vector valued function F. Setting F = Vu gives what 
is called Green's First Identity, 

/ / (GX7 2 u + Vm • VG) dA = [ G(Vu-h)dS (2) 
J Jd Jc 

Interchanging G and u and subtracting gives Green's Second Identity, 

/ / (uV 2 G - GV 2 u) dA = [ (uVG-GVu)-hdS. (3) 
J Jd Jc 

2 Solution of Laplace and Poisson equation 

Ref: Guenther & Lee, §5.3, §8.3, Myint-U & Debnath §10.2 - 10.4 
Consider the BVP 

V 2 m = F in D, (4) 
u — f on C. 

Let (x, y) be a fixed arbitrary point in a 2D domain D and let (£, if) be a variable 
point used for integration. Let r be the distance from (x,y) to (£,rj), 



r = \I — x) 2 + (r) — y) 2 . 



Considering the Green's identities above motivates us to write 

V 2 G = 6(Z-x,r t -y) = 6(r) in D, (5) 
G = on C. 



The notation 5 (r) is short for 5 (£ — x,r) — y). Substituting (4) and (5) into Green's 
second identity (3) gives 

u(x,y)- [ [ GFdA = [ fVG ■ MS 
J Jd Jc 

Rearranging gives 

u(x,y)= [ [ GFdA + I fVG ■ hdS (6) 
J Jd Jc 

Therefore, if we can find a G that satisfies (5), we can use (6) to find the solution 
u (x, y) of the BVP (4). The advantage is that finding the Green's function G depends 
only on the area D and curve C, not on F and /. 

Note: this method can be generalized to 3D domains. 



2.1 Finding the Green's function 

To find the Green's function for a 2D domain D, we first find the simplest function 
that satisfies V 2 f = 5 (r). Suppose that v (x,y) is axis-symmetric, that is, v — v (r). 
Then 

V 2 v = v rr H — v r = S (r) 
r 

For r > 0, 

v rr H — v r = 
r 

Integrating gives 

v = A\nr + B 

For simplicity, we set B = 0. To find A, we integrate over a disc of radius e centered 
at (x,y), D e , 

1=1 [ 6(r)dA= f [ V 2 vdA 




From the Divergence Theorem, we have 

V 2 vdA= I Vv-ndS 

D £ Jc £ 

where C £ is the boundary of D e , i.e. a circle of circumference 1-ne. Combining the 
previous two equations gives 

1=1 Vv-ndS = [ ^ dS= f -dS = 2rrA 
Jc, Jc, or r=p Jn e 



<C e JC, 

Hence 



v (r) = — lnr 



This is called the fundamental solution for the Green's function of the Laplacian on 
2D domains. For 3D domains, the fundamental solution for the Green's function of 
the Laplacian is — l/(4irr), where r = \J (x — £) 2 + (y — r/) 2 + (z — () 2 . 

The Green's function for the Laplacian on 2D domains is defined in terms of the 
corresponding fundamental solution, 

G(x,y;C,v) = ^~ lnr + /i ' 
h is regular, 

V 2 h = 0, (^)GA 

G = (^)eC. 

The term "regular" means that h is twice continuously differentiate in (£,77) on D. 
Finding the Green's function G is reduced to finding a C 2 function h on D that 
satisfies 

V 2 h = (^)eD, 

h = -i-lnr ((.^C. 

The definition of G in terms of ft, gives the BVP (5) for G. Thus, for 2D regions D, 
finding the Green's function for the Laplacian reduces to finding h. 

2.2 Examples 

Ref: Myint-U & Debnath §10.6 

(i) Full plane D = R 2 . There are no boundaries so h = will do, and 

G=^lnr=i-ln [(£ - x) 2 + (r, - y) 2 } 

(ii) Half plane D = {(x,y) : y > 0}. We find G by introducing what is called an 
"image point" (x,—y) corresponding to (x,y). Let r be the distance from (£,77) to 
(x,y) and r' the distance from (^,77) to the image point (x, —y), 



r = y - x) + (77 - y) , r' = ^ - x) + (77 + y)' 

We add 

1 , , 1 



h = ~y lnr' = ~y In y (f - x) 2 + (77 + y) 2 

to G to make G = on the boundary. Since the image point (x, —y) is NOT in D, 
then /i is regular for all points (£,77) € -D, and satisfies Laplace's equation, 

V 2 h _^ —-0 
<9£ 2 (9r/ 2 
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for (£,77) G -D. Writing things out fully, we have 

G = ^lnr + / i =^lnr-^lnr' = ^ln^ = ^ln^4±^4 (7) 

G (x, y; £, 77) is plotted in the upper half plane in Figure 1 for (x, y) = (y/2, y/2) . 
Note that G — > —00 as (£,77) — > (x,y). Also, notice that G < everywhere and 
G = on the boundary 77 = 0. These are, in fact, general properties of the Green's 
function. The Green's function G (x, y; £, 77) acts like a weighting function for (x, y) 
and neighboring points in the plane. The solution u at (x, y) involves integrals of 
the weighting G (x, y; £, 77) times the boundary condition / (£, 77) and forcing function 

On the boundary C, 77 = 0, so that G = and 



VG • n = - — 

orj 



y 



The solution of the BVP (6) with F = on the upper half plane D can now be 
written as, from (6), 

u (x, y)= f fVG -MS= V - r —L®—dS, 



which is the same as we found from the Fourier Transform, on page 13 of fourtran.pdf. 
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(iii) Upper right quarter plane D = {(x,y) : x > 0,y > 0}. We use the image 
points (x,-y), (-x,y) and (-x,-y), 



G = 



1 

2tt 



In \J(£- xf + (77 - yf - — In ^/(f - a;) 2 + (77 + y) 



(8) 



For (£, 77) G C = dD (the boundary), either £ = or 77 = 0, and in either case, G — 0. 
Thus G = on the boundary of -D. Also, the second, third and fourth terms on the 
r.h.s. are regular for (£, 77) G -D, and hence the Laplacian V 2 = d 2 /d^ 2 +d 2 /drj 2 of each 
of these terms is zero. The Laplacian of the first term is 5 (r). Hence V 2 G = 5(r). 
Thus (8) is the Green's function in the upper half plane D. 
For (£,77) £ C = dD (the boundary), 



fVG ■ MS 



a 



00 
00 



7(0,77) 



dG 



dG 



^+ / /(£,o) 

#77 



dG 
drj 



d£ 



ij=0. 



dr) 



r]=0 



Note that 



dG 

dG 
drj 



4yxr] 



T]=0 



71 (x 2 + (y + r]) ) (x 2 + (y - 77) ) 
^((^-0 2 + y 2 ) {{x + if + y 2 ) 



The solution of the BVP (6) with F = on the upper right quarter plane D and 
boundary condition u = f can now be written as, from (6), 



u(x,y) 



fVG ■ MS 



c 
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vf (0,77) 
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a (CO) 
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drj 
di 



(iv) Unit disc D = {(x,y) : x 2 + y 2 < 1}. By some simple geometry, for each 
point (x,y) G D, choosing the image point (x',y') along the same ray as (x,y) and a 
distance 1/ ^ x 2 + y 2 away from the origin guarantees that r/r' is constant along the 
circumference of the circle, where 



(£-x) 2 + (77 -yf, 



(Z-xf + in-yf. 
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[DRAW] Using the law of cosines, we obtain 



p 2 + p 2 - 2ppcos (o - 0) 

('-•) 



r' 2 = p 2 + 2- cos 
P P 

where p = \/ x 1 + y 2 , p = a/£ 2 + rf and 0, 6 are the angles the rays (x, y) and (£, rf) 
make with the horizontal. Note that for (£, 77) on the circumference (£ 2 +^ 2 = p 2 = 1), 
we have 

r 2 1 + p 2 -2pcos (fl- 9 



P\ 



P 



l + ^-2jcos^-^ 
Thus the Green's function for the Laplacian on the 2D disc is 

I p 2 + p 2 - 2ppcos (e- 9] 
In 



G(£,ri;x,y) 



1 r 
In — 



Note that 



VG-n 



2tt r'p 4tt p2;5 2 + ! _ 2p p C0S ( - 9 

1-p 2 



<9G 

<9p 



1 



2tt 



p=i ~" 1 + p 2 -2pcos 
Thus, the solution to the BVP (5) on the unit circle is (in polar coordinates), 

r-27T -1 9 



U(p,0) 



1 

27 



1-P 2 



dO 



1 + p 2 - 2p cos [e - 0) 

p 2 + p 2 - 2pp cos ( - 6> 



2 p2n pi 

4tt Jo io 



+ 



In 



p 2 p 2 + 1 - 2ppcos 0-0 



[p, 9 s ) pdpdO 



The solution to Laplace's equation is found be setting F — 0, 



u(p,0) = 



2tt 



2tt 



1 



l + p 2 -2pcos 0-0 



dO 



This is called the Poisson integral formula for the unit disk. 



2.3 Conformal mapping and the Green's function 

Conformal mapping allows us to extend the number of 2D regions for which Green's 
functions of the Laplacian V 2 w can be found. We use complex notation, and let 
a = x + iy be a fixed point in D and let z = £ + irj be a variable point in D (what 
we're integrating over). If D is simply connected (a definition from complex analysis), 
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then by the Riemann Mapping Theorem, there is a conformal map w (z) (analytic 
and one-to-one) from D into the unit disk, which maps a to the origin, w (a) = 
and the boundary of D to the unit circle, \w (z)\ — 1 for z G dD and < \w (z)\ < 1 
for z G D/dD. The Greens function G is then given by 



G 



2tt 



In In; (z) 



To see this, we need a few results from complex analysis. First, note that for z G dD, 
\w (z)\ = so that G — 0. Also, since w (z) is 1-1, \w (z)\ > for z ^ a. Thus, we 
can write u> (z) = (z — a) n H (z) where H (z) is analytic and nonzero in D. Since 
w (z) is 1-1, \w' (z)\ > on D. Thus n — 1. Hence 



iu (z) — (z — a) H (z) 



and 



where 



|z - 
1 

2^ 



2vr 



lnr + /i 



\n\H(z) 



Since H (z) is analytic and nonzero in D, then (l/27r) lnif (z) is analytic in D and 
hence its real part is harmonic, i.e. h = 3? ((l/27r) In H (z)) satisfies V 2 h = in D. 
Thus by our definition above, G is the Green's function of the Laplacian on D. 
Example 1. The half plane D = {(x,y) : y > 0}. The analytic function 



w (z) 



a 



or 



maps the upper half plane D onto the unit disc, where asterisks denote the complex 
conjugate. Note that w (a) = and along the boundary of D, z = x, which is 
equidistant from a and a*, so that \w (z)\ = 1. Points in the upper half plane (y > 0) 
are closer to a = x + iy, also in the upper half plane, than to a* = x — iy, in the 
lower half plane. Thus for z G D/dD, \w (z)\ = \z — a\ / \z — a*\ < 1. The Green's 
function is 

G — — In |k; (z) I = — In r = — In — 



2tt 1 v " 2tt \z 
which is the same as we derived before, Eq. (7). 
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3 Solution to other equations by Green's function 

Ref: Myint-U & Debnath §10.5 

The method of Green's functions can be used to solve other equations, in 2D and 
3D. For instance, for a 2D region D, the problem 

V 2 u + u = F in D, 
u = f on dD, 

has the fundamental solution 




where Y (r) is the Bessel function of order zero of the second kind. The problem 

V 2 m-m = F in D, 
u = f on dD, 

has fundamental solution 

where K (r) is the modified Bessel function of order zero of the second kind. 

The Green's function method can also be used to solve time-dependent problems, 
such as the Wave Equation and the Heat Equation. 
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